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Several authors, most notably Heckler, have claimed that the observable Hawking emission from 
a microscopic black hole is significantly modified by the formation of a photosphere around the 
black hole due to QED or QCD interactions between the emitted particles. In this paper we 
analyze these claims and identify a number of physical and geometrical effects which invalidate 
these scenarios. We point out two key problems. First, the interacting particles must be causally 
connected to interact, and this condition is satisfied by only a small fraction of the emitted particles 
close to the black hole. Second, a scattered particle requires a distance ~ Ejm\ for completing 
each bremsstrahlung interaction, with the consequence that it is improbable for there to be more 
than one complete bremsstrahlung interaction per particle near the black hole. These two effects 
have not been included in previous analyses. We conclude that the emitted particles do not interact 
\f*l • sufficiently to form a QED photosphere. Similar arguments apply in the QCD case and prevent a 

■ QCD photosphere (chromosphere) developing when the black hole temperature is much greater than 
Aqcd, the threshold for QCD particle emission. Additional QCD phenomenological arguments rule 
out the development of a chromosphere around black hole temperatures of order Aqcd- In all cases, 
the observational signatures of a cosmic or Galactic halo background of primordial black holes or 
an individual black hole remain essentially those of the standard Hawking model, with little change 

q , to the detection probability. We also consider the possibility, as proposed by Belyanin et al. and 

■ D. Cline et al., that plasma interactions between the emitted particles form a photosphere, and we 
' conclude that this scenario too is not supported. 
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In 1974 Hawking showed that black holes should continually emit radiation as a natural consequence of quantum field 
theory in curved spacetime P, Q ■ Hawking radiation is analogous to thermal radiation from a black body and leads 
Q\ • to the identification of a classical black hole temperature. In the case of a Schwarzschild black hole, this temperature 
is inversely proportional to the black hole mass and therefore increases as the black hole radiates [l| • Only black holes 
with mass below 10 26 g have a temperature today greater than that of the cosmic microwave background. Such black 
holes may have been produced in the early Universe with dramatic cosmological consequences. Formation scenarios for 
primordial black holes (PBHs) include the collapse of overdense regions due to primordial inhomogeneities, especially 
J^j ' those generated by inflation, a softening of the equation of state or bubble collisions at cosmological phase transitions, 
and the collapse of oscillating cosmic string loops or domain walls. The cosmological consequences include modifying 

■ the usual picture of cosmological nucleosynthesis, distorting the spectrum of the microwave background radiation, 
generating a cosmic baryon asymmetry, removing magnetic monopoles, and contributing to the dark matter. Even 
if PBHs never existed, studying them therefore places important constraints on models of the early Universe. All of 
these formation scenarios and cosmological consequences have recently been reviewed by Carr 

Of particular interest are those PBHs with initial mass ~ 10 15 g which are evaporating at the present epoch, since 
these should contribute to the observable 7-ray, cosmic ray and neutrino backgrounds (Q and references therein). 
Indeed the detection of such backgrounds or of individual PBH bursts would be the first confirmation of the funda- 
mental relationship linking black holes, general relativity and thermodynamics (see also Ref. To date, no such 
PBH burst or background has been conclusively detected, although D. Cline et al. [|| @, Q have attributed some 
short-timescale gamma-ray bursts to PBH explosions. However, observations of the 100 MeV gamma-ray background 
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imply that 10 15 g PBHs must have Cl p bh ^ 10 -8 where fl p bh is their present density in units of the critical density. 
The 5-year data from the Whipple gamma-ray telescope impose the strongest current limit on the PBH explosion 
rate . With the Standard Model of particle physics, this is much weaker than the gamma-ray background limit on 
flpbh if the PBHs cluster smoothly within the Galactic halo, but it could be stronger if PBHs were highly clumped. 
The £lpbh limits from the Galactic cosmic rays are comparable with the gamma-ray background limit but depend on 
both the PBH clustering and the cosmic ray propagation model [ToL [Tlj|. 

Black holes of ~ 10 15 g have an initial temperature of about 20 MeV. Because a black hole will emit all particle 
species with rest mass less than or of order the black hole temperature, the emission from a 20 MeV black hole 
is sensitive to the threshold for quark and gluon production, Aqcd ~ 200 — 300 MeV. In the standard Hawking 
model [l2| , the particles directly evaporated from a black hole are those which appear noncomposite compared to the 
wavelength of the radiated energy (or equivalently the black hole size) at a given temperature. These particles can then 
decay or form composite particles as the emission streams away from the black hole. A black hole with a temperature 
somewhat below Aqcd should directly emit neutral and charged pions, as well as photons, electrons, positrons, muons, 
neutrinos and gravitons. Once the energy distribution of the Hawking emission significantly extends above Aqcd, 
the black hole should directly emit quarks and gluons rather than pions, with the quarks and gluons fragmenting 
and hadronizing in jets into the astrophysically stable species (7, p, p, e + , e~ , v and 77) after emission. This process 
is analogous to QCD jet decay in accelerators. By convolving the Hawking evaporation with the HERWIG Monte 
Carlo jet code, MacGibbon and Webber found that the particle flux from 0.3 — 100 GeV black holes is predominantly 
comprised of jet fragmentation products [T2I ]. They found that, after decay, the flux peaks at an energy of around 
100 MeV, roughly independent of the black hole temperature, with the average energy of the final species scaling 
more weakly than the black hole temperature. In contrast, prior to decay, the peak and average energies of the direct 
Hawking emission arc proportional to the black hole temperature. 

In the Heckler photosphere scenario, as originally described in Refs. [H, [l4j|, the situation is radically different. 
Above a certain black hole temperature, the particles are alleged to interact sufficiently after emission to form a 
quasithermal photosphere at a distance a -4 times the Schwarzschild radius from the black hole [l3j], where a = 
e 2 /he 1/137 is the fine structure constant. In this model, the relevant QED interactions are bremsstrahlung and 
electron-photon pair production. The critical temperature above which this happens is estimated to be T cr u ~ 45 

5 /2 

GeV. When QCD effects are included, a lower transition temperature T cr u <~ Aqcd/o-s is derived where a s is 
the strong coupling constant: a s m 0.12 at energy scales around 100 GeV but it increases at lower energies and its 
precise behavior is unknown around Aqcd- Once the photosphere develops, the high energy emission is reprocessed 
to lower energies, drastically reducing the flux of high energy particles and significantly weakening the detectability 
of individual high-temperature black holes. Because the best limits on Vt p bh are derived from the present 100 MeV 
emission from a Galactic or extragalactic background of PBHs with temperatures well below T cr n, the limits on 
£l p bh are weakened only slightly by photosphere development [l4j|. However, photosphere development would have 
considerable implications for experimental efforts to detect ultra-high-energy cosmic rays from the final explosive 
phases of PBHs in the present epoch [l^]. Indeed, this has led to decreased motivation among groups searching for 
such signatures. 

J. Cline, Mostoslavsky, and Servant [l6| have numerically explored the development of the photosphere in the 
original Heckler model. They conclude that photosphere formation should set in above a slightly lower value of T cr n, 
weakening the bound on Q p bh still further and significantly changing the burst spectra. These authors, however, 
assume the same initial description of interactions as Heckler [l3l . [l4| , so their work is not an independent derivation 
or confirmation of the model. Kapusta [l?], [Ul| an d Daghigh and Kapusta [ljj [2^, [2l| have also performed calculations 
related to Heckler's model. 

In this paper we will raise a number of problems with the Heckler scenario. These include various geometrical 
effects, the form of the plasma mass correction, the validity of the perfect fluid assumption and QCD phenomenological 
considerations. However, the two most serious omissions in the Heckler scenario are as follows: (i) two particles must 
be in causal contact to interact, and this is satisfied only by a small fraction of the emitted rclativistic particles 
when they are sufficiently close to the black hole that they would otherwise have significant interactions; and (ii) in a 
bremsstrahlung interaction the electron initially forms off-shell and must travel a distance in the black hole frame much 
greater than ~ l/m e before completing each interaction and becoming on-shell. If it is 'bumped' by encounters with 
other particles within ~ l/m e of the hole, the Landau-Pomeranchuk-Migdal (LPM) effect implies that the particle 
generally undergoes at most one completed bremsstrahlung interaction as it streams away from the black hole. Thus 
we conclude that a persistent QED or QCD photosphere cannot form around a black hole. Rare interactions could 
slightly modify the 7, p, p, e + , e~ , v and v signatures of PBHs but not significantly. 

Other authors, such Belyanin, Kocharovsky, and Kocharovsky [23] and D. Clinc ct al. @, have proposed 
photosphere models which are unrelated to the Heckler mechanism. These models arc mainly motivated by an 
attempt to explain certain short-period gamma-ray bursts. After analyzing the Bclyanin-Kocharovsky-Kocharovsky 
and D. Cline et al. models, we conclude that particle interactions are also insufficient to form a QED or QCD 
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photosphere in these scenarios. However, the relevance of PBH explosions at the present epoch should not be ruled 
out as some other PBH mechanism which has not yet been considered might explain the observed bursts. 

The plan of this paper is as follows. In Section II we review the fundamentals of Hawking evaporation. In Section 
III we outline the original Heckler photosphere scenario. In Section IV we analyze this scenario in detail. In Section V 
we summarize the results of our analysis and their observational implications. Other photosphere models are discussed 
in Section VI. Our conclusions are presented in Section VII. Except in Sections II and VIA, we will use Planck units 
in which h = k = G = c = Aits® = 1 . 



II. HAWKING EVAPORATION 



A Kerr-Newman black hole with angular velocity Q, and electric potential <& radiates particles with total energy 
between (E, E + dE) at a rate [J, [3| 



dN s = Y, TsnldE 
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(2.1) 



per degree of particle freedom. Here s is the particle spin, nh is its axial quantum number or angular momentum, I 
is its orbital angular momentum, q is its charge and T sn i is the dimcnsionlcss absorption probability for the emitted 
species. In general, T sn i is a function of E, fl, $ and the surface gravity of the black hole k, as well as the internal 
degrees of freedom and rest mass of the emitted particle. Equation (|2.1|) implies that a charged, rotating black 
hole will preferentially emit particles which have the same sign of charge and spin as the black hole. The angular 
momentum of a black hole is emitted several times faster than its mass I23|. An electrically charged black hole quickly 
discharges provided M bbj < 10 5 M Q , where M bb is the black hole mass [24l I25I I2II I27L l28 |. Thus we will henceforth 
assume an uncharged, nonrotating Schwarzschild black hole, i.e. f2 = $ = and n = c 4 / (4GM b h)- 
The temperature of an uncharged, nonrotating black hole is [l[ 

kT bh = Jgj- = 1.06 (^-Y 1 GeV. (2.2) 



8irGM bh ■ Vl° 13 £ 

At all temperatures, a black hole will radiate the massless elementary particles: s = 1/2 neutrinos (if we ignore 
neutrino mass), s = 1 photons and s = 2 gravitons. Massive elementary particles are evaporated in significant 
numbers once the peak in the energy distribution given by Eq. (|2.ip is of order their particle rest mass. Since the 
exponentially damped tail in Eq. (|2.ip extends to infinite energies, a nonzero contribution from any massive species 
is always present to some extent. 

At high energies, the sum over n and I of the absorption probabilities for both massless and massive species 
approaches the geometric optics limit [H, [53, [H| : 

T s (M bh ,E) ^Y, T ^ * 2?G ^p 2 . (2.3) 
At low energies, T s is suppressed for massless s = and s = 1/2 species: 

1 s =0 {M bh , E ) W , 1 a _jy 2 {M bh , E ) W ^6 ■ \ ZA ) 



For s = 1 and s = 2 massless bosons, the E —> fall-off is steeper: 

64:G 4 M* h E 4 „ _ 256G 6 M^E 6 

3fr 4 c 



l s=1 {M bh ,E) -— -35 , V s=2 {M bhl E) 7^6~i8 ■ V 2 - 5 ) 



The absorption probability for nonrelativistic massive particles is more strongly suppressed at low energies than 
for massless or relativistic particles. For (Mbft,/10 13 g) (///GeV) < 10, where fi is t he p article rest mass energy, 
T s=1 / 2 (M b h, E,fi) at E = /i differs from the corresponding [i, = value by up to 45% [28|]. More importantly, since 
E is the total energy, the spectrum of a massive species must be zero for E < \i. These effects have a strong damping 
influence on the emission spectrum when the peak energy is close to the particle rest mass threshold. 

Prior to particle decays, the peaks in the flux per particle mode measured at infinity occur at [23|, [2?], HH, 



E s=0 » 2.81T 6h , E s=1/2 = 4.02T bh , E s=1 = 5.77T bh (2.6) 
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and the total instantaneous flux emitted by a black hole is 

dN v-^ f°° dNi , , , 

where n^ is the number of degrees of freedom (or modes) per particle species i. The standard Glashow-Weinberg-Salam 
model with three generations has 2 modes per s = 1/2 neutrino or antineutrino, 4 modes per s = 1/2 e^, /i^ or 
lepton. 12 modes per s = 1/2 quark flavor, 16 modes per s = 1 gluon, 2 modes per s = 1 photon, and 2 modes per 
s = 2 graviton. With 40 — 100 elementary modes for T^h ~ 0.3 — 100 GeV, the total emission rate of particles (prior 
to their decays) is 

— ~ 10 ' 2c3 (2 &) 

dt ~ GM bh ' { ' ' 

Combining Eq. (|2.7|) with the particle energy gives the instantaneous mass loss of the black hole. This can be 
integrated to derive the lifetime of a black hole of initial mass Mj (30j : 

r evap = 6.24 x 10~ 27 Mff (Mi)' 1 s. (2.9) 

The factor / (Mi) depends on the number of emitted species and is normalized to unity when only massless species 
are emitted. Following the method of Ref . [3(3| , we can now update the value of M» , the initial mass of a black hole 
whose lifetime is the present age of the Universe. Using the 3-year Wilinson Microwave Anisotropy Probe (WMAP) 
best-fit value for the age of the Universe in the Lambda Cold Dark Matter (LCDM) model, t u = 13.73+q 'j^ x 10 9 yr 
[3ll | . we find that 

M* = (5.00 ± 0.04) x 10 14 g. (2.10) 

To aid our conceptual discussion of interactions between the emitted particles, we note the magnitudes of various 
physical quantities. The radius of the black hole is 

2GMbh -, , n 1ri -l5 ( \ in-15 ( Tbh \ ;„„\ 

r bh = ~^T- = 1-49 x 10 jcm= 1.57 x 10 ^— J cm. (2.11) 

For directly emitted relativistic particles of energy E, the reduced de Broglie wavelength (a measure of the effective 
interaction range for interactions which fall off as some power of 1 / E) is 

A 5 = J = 0.197 x nr 13 (tJ^) cm. (2.12) 

In the nonrelativistic case, the factor E/c in Eq. I|2.12[) is replaced by the momentum. Since the flux per s = 0, 1/2 
or 1 mode peaks at E pea k ~ (3 — 6) T^u, we have 

Aft = (3 - 7) x 10- 15 (Jr^) cm« (2-5)rw». (2.13) 

A typical emitted particle therefore has a wavelength comparable to or somewhat larger than the size of the hole 2rft^. 
Additionally, the particles can be considered 'evaporated' at a distance r ev ks Srbh/Z since the outgo ing solution of 
the Hawking wave equation here already closely approximates the solution for the particle at infinity [32j. 

If the time between successive emissions from the black hole is less than Aft /c, or equivalently if the distance travelled 
by a particle before the hole emits another particle is less than Aft, we would expect interactions between the particles 
regardless of their species. Equation (|2.8[) gives the average time between successive emissions of elementary particles 
as 

lOOGMft^ 

At« ^— — , (2.14) 



so 



c Ai«J^»Aft (2.15) 
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and such interactions are unlikely. More precisely, Olicnsis and Hill have calculated numerically that more than 99% 
of the particles emitted over a hole's lifetime satisfy the condition At > h/E [33|. This confirms the appropriateness 
of using E peak in Eq. (|2.15|> . 

Similarly, from Eq. (|2.2p . once the temperature is above a few times Aqqd, the energy scale on which the QCD 
force becomes strong, the time between successive emissions is much shorter than the QCD timescale H/Aqqd, and 
the size of the black hole 2rbh is much smaller than the associated distance Hc/Aqcd- Thus, for Tbh >> Aqcd, the 
Hawking evaporation of a colored particle is unaffected by the emission of other colored particles and the particle can 
be regarded as asymptotically free. On the scales associated with Aqcd in the black hole rest frame, the emission of 
successive colored particles at these temperatures appears as effectively simultaneous production of highly energetic 
(E » Aqcd) quarks and gluons at a point. This is analogous to QCD jet production in e + e~ annihilation accelerator 
events. As the quarks and gluons propagate away from the black hole, they should fragment into further quarks and 
gluons and finally cluster into colorless hadrons once they have travelled a distance (Hc/Aqcd) [E a / Ha) from the black 
hole. Here E a represents an 'average' energy of the jet particles involved in the decay and fi a is an average 'effective' 
rest mass energy. 

At temperatures above Aqcd, the observable spectra from the black hole, found by convolving the Hawking 
formula with jet decays, are dominated by the QCD decay products. Once the emitted particles have decayed into 
the astrophysical stable species, the average energy for Tbh ~ 0.3 — 100 GeV scales as approximately T^ 2 for 7, e + , 
e~, v and 77, and X^ 8 for p and p, rather than as Ty, as given by Eq. (|2.6p [l2| . The particle flux after decay, dN/dt, 
for ~ 0.3 — 100 GeV scales as approximately T 1 ^ 6 [12||, rather than as Tbh as implied by Eq. (|2.8p . 



III. THE HECKLER SCENARIO 

We first focus on the QED interactions in the Heckler scenario, as originally derived in Refs. 0, E3- Consider a 
Tbh >> m e black hole which is emitting Hawking radiation in random radial directions. From Eq. (|2.ip the number 
density at radius r of electrons and positrons which were directly radiated by the black hole (and not the result of 
particle decays) is 

10" 4 

where we are now using units with h=k = G = c= \. If the number density is sufficiently high, the electromagnetic 
particles could lose energy via two-body bremsstrahlung e + e^e + e + 7 and pair production e + 7 — * e + e + + e~ 
as they propagate away from the hole. Each such interaction would increase the number of electromagnetic particles 
by a factor of 3/2. 

To calculate the number of scatterings per particle, Ref. \\M ta kes the bremsstrahlung cross-section in the center- 
of-mass frame — averaged over the photon energy us — to be |33l . l34j 

1 f E d<j brem 8a 3 2E 

CTbrem = -=; / U) : duj « ^-lll . (3.2) 

E J dui ra e m e 

The center-of-mass frame for each interaction is assumed to be the rest frame of the black hole. In Eq. (|3.2p . E ss 4Tbh 
is the energy of the initial electron or positron. Since the average momentum exchanged by the virtual photon in the 
bremsstrahlung interaction is ~ m e , the Heisenberg Uncertainty Principle implies that the two initial particles must 
be within ~ l/m e of each other to interact significantly. In Ref. [l3| the 'formation' distance for the final on-shell 
electron and photon is also assumed to be ~ l/m e in the center-of-mass frame of the interaction. In the relativistic 
limit, the bremsstrahlung cross-sections for two incident electrons or two incident positrons or an electron and a 
positron are the same. The pair-production cross-section can be derived from the bremsstrahlung cross-section using 
the substitution law (see Section IV. K) and has the same functional form in the relativistic limit [351 ] - 

In Ref. [IH , m e in Eq. (|3.2[) is replaced by an effective electron mass m' e which is taken to be the vacuum electron 
mass m e augmented by a 'plasma mass' m pm to account for Coulomb screening by the background of other nearby 
emitted electromagnetic particles. The electron mass is approximated by the correction applicable for an isotropic 
thermal plasma 



;— -— — 2 _ Anan (r) 

e ' pmi m pm ~ p i 

which includes the further approximations (1/E) f» l/E av w 1/ATbh- This increases the effective electron mass, 
resulting in a decreased cross-section, within a distance ~ O.la 1 / 2 ^." 1 of the black hole. However, there is an 
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inaccuracy here in that the emitted particles do not form an isotropic thermal plasma, as pointed out by Heckler [131 ] 
and Cline, Mostoslavsky, and Servant [l(| themselves. 

The number of scatterings, M , experienced by a particle by the time it reaches radius R is found by calculating 

with A (r), the mean-free-path of an individual particle, defined as 

A (r) = (n (r) ov e m iVei) -1 (3.5) 

and 

2 J no(r). (3.6) 

The factor of 3/2 accounts for the increase in the number of electromagnetic particles at each bremsstrahlung interac- 
tion, and the relative velocity between the interacting particles is taken to be c. Only the lower limit in integral (|3.4|) 
is important. As r — > oo, the number density of the emission decreases as r~ 2 and the scattering cuts off. Heckler first 
neglects the plasma mass correction and takes the lower limit in integral (|3.4j) to be l/m e . In this case, the integral 
becomes 

N ~ In , (3.7) 

27r 4 m e m e 

which exceeds 1 for 

Tbh > -^m e ~ 20 TeV. (3.8) 

cr 

The analysis is then extended to include the expression for the plasma mass given by Eq. (|3.3[) and the lower limit 
in integral (|3.4p is taken to be the radius of the black hole. For r small enough that m pm (r) > m e , the number of 
scatterings is damped by the plasma mass. The dominant contribution to the integral then comes from the value of 
r at which m pm ps m e and this gives 

^«T7J — ln — . ( 3 - 9 ) 

which exceeds 1 for 

T ^i^ji K45GCV - (3 - 10) 

On numerically integrating Eq. (|3.4[) . Heckler finds that AT = 1 when T cr a « 45.2 GeV. Thus, although the cross- 
section abrem ot a 3 remains small, the growth in number density around the black hole at high T^u leads to significant 
interaction in this model. 

The region for which TV >> 1 is the photosphere region. It is claimed that when >> T cr n the multiple 
scatterings are sufficient to produce an interacting quasithcrmal fluid, which flows away from the black hole. The 
photosphere is found to have inner and outer radii given by 

An 16n2 ~m9 ~ 1 f T "hY /2 

-r bh ps 10 J r bh , r out ps — — — . (3.11) 



Approximating this region by a perfect fluid, the average energy of the particles at the outer boundary of the 
photosphere, in effect the observable average energy of the emitted particles, is estimated to be 

E bs ~ m e {T bh /T crlt ) 1/2 , (3-12) 

which is considerably less than the energy E av ps 4T b h expected without photosphere development. (These values do 
not include the electrons and photons produced by the decays of other particles radiated by the black hole.) Therefore 
the flux of photons observed from a black hole with a QED photosphere is 



dN 1 ^ T b h ( Tbh x 



dt a 2 IT, 



(3.13) 
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neglecting particle decays, compared to dN 1 /dt ps 10~ 2 Tbh without a photosphere. 

Extending the Heckler model to QCD effects, Eq. (|3.2p is replaced in Ref. [l3[ by the cross-section for the gluon- 
brcmsstrahlung process q + q^q + q + g: 

O-ferem « —fin (3.14) 

771^ TO, 

where a s is now the strong coupling constant. The effective quark mass, m q , is also enhanced by a plasma term: 



< = Alan + ^ <{MQCD i ~ l? ' • (3-15) 

1 QCD ( 47r 2 r )2 



Proceeding as before, one arrives at a rough estimate for the onset of a photosphere due to QCD interactions: 
Nqcd >> 1 when T^h exceeds 

T? r T *A Q c D /a^>A QCD . (3.16) 

The precise value of T® rit D is sensitive to threshold effects and the energy dependence of a s around Aqcd- Because 
Tbh must be greater than about Aqcd for a black hole to emit QCD particles, T® rit D cannot be less than about 
Aqcd- 

Cline, Mostoslavsky, and Servant [lj| adapted a code developed to study heavy ion collisions to obtain numerical so- 
lutions of the Boltzmann equation in the context of Heckler's model. As discussed in Section IV. L, these authors found 
that QED and QCD photospheres may form but disagreed with some features of the original Heckler model. Once the 
photosphere formed, they also considered the effect of elastic Compton scattering (which conserves particle number) 
and found that this can extend the photosphere when > 5 TcV. However, because Cline, Mostoslavsky, and 
Servant employed the same initial description of interactions as Heckler, their work is not an independent derivation 
or confirmation of photosphere development. 



IV. IS THE HECKLER SCENARIO CORRECT? 



In this section, we show that several aspects of the original Heckler model arc invalid and conclude that QED 
photosphere development is ruled out. We begin by focusing on the QED brcmsstrahlung interactions. Extension 
to QED pair production, as well as the corresponding QCD interactions, is straightforward and discussed in Section 
IV. K. As in the original Heckler model, we will consider all emitted particles to be initially moving in random radial 
directions away from the black hole. The effect of black hole size on the emission directions is negligible for r » r^h- 
We will also ignore general relativistic effects, which are negligible for r > lOr^. 



A. Ubrem in More Detail 



First we consider whether the bremsstrahlung cross-section used by Heckler is appropriate. The double-differential 
cross-section for two-body electron bremsstrahlung is given in Ref. [HJ . In this interaction, two electrons and/or 
positrons exchange a virtual photon, forcing one electron off-shell. The off-shell electron then decays into an on-shell 
electron by emitting a photon. The eight relevant Feynman diagrams for lowest order perturbation theory are shown 
in Fig. [1] .'Hi . The cross-section Eq. (|3.2[) applies in the relativistic limit for e~e~, e + e + and e~e + bremsstrahlung 
because the interactions associated with the exchange diagrams, Figs. 1 (e), (f), (g) and (h), are suppressed at high 
energy. (In the nonrelativistic limit, the e~e + bremsstrahlung cross-section is greater than the e~e~ and e + e + cross- 
sections due to the dipole nature of e~e + [H[.) In the relativistic limit, the average angle between the final on-shell 
electron and photon is [3J] 

l*L - 1 ( 4 -i) 

in the center-of-mass frame of the interaction, where E is the energy of each initial electron. The average final energies 
of the on-shell electron and photon are [34| 

E e « lu w £72. (4.2) 
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FIG. 1: The eight Feynman diagrams for electron-electron bremsstrahlung in lowest order perturbation theory. The labels, 
which are defined in the text, correspond to the 4-momenta of the particles involved in the interaction. 



Because <p av is so small, the majority of interactions involve very small momentum transfer, typically of order m e . 
Thus the directions of the electrons are modified little by the interaction but their energy loss can be substantial. 

Since the average energy of the final photon is of order E, it is valid to use the energy-averaged two- body 
bremsstrahlung cross-section given by Eq. (|3.2[) which subsumes the infrared divergence as ui — > 0. In Ref. [3_4j, 
d(Jbrem / dtu is evaluated by numerically integrating the complicated double-differential cross-section. The cross-section 
does not contain the collincar divergences which would arise at <fi = 0° and <f> = 180° if factors of order m e were not 
correctly included in the derivation. In the non-relativistic limit, abrem decreases more strongly with decreasing E 
than Eq. (|3.2[) Q and becomes process-dependent [HjHJ, as noted above. 

It should be stressed that the Heckler effect is not due to an infrared or collincar divergence in the cross-section 
but to the high number density of emitted particles flowing isotropically and relativistically away from the black hole. 
As a consequence of the emission distribution, the center-of-mass frames for most of the interactions are not moving 
highly relativistically with respect to the black hole, as we next show. 



B. Non-interaction of Radially Comoving Particles 



In Section II, we saw that one should not expect QED or QCD interactions between emitted particles if the 
interaction scale falls off as 1/E (i.e. if the cross-section falls off as 1/E 2 ). In the Heckler scenario a quasithermal 
photosphere develops from QED and QCD bremsstrahlung and pair-production interactions, which do not obey this 
constraint. We now demonstrate that particles moving in similar radial directions do not interact significantly even 
via bremsstrahlung or pair production. 

Consider a typical directly emitted particle of rest mass /j, << Tbh and energy E sa 4Tbh at a distance r >> rbh 
from the black hole. At this distance, special but not general relativistic effects are relevant. In the rest frame of the 
black hole, the particle has an individual Lorentz factor 



7 « {2-KMbhfi) 



(4.3) 



FIG. 2: Illustrating how a particle emitted by a black hole of radius rth can interact with other particles emitted in a similar 
radial direction where the angle Oth is the deviation from the radial direction due to the finite size of the black hole. 



and velocity 

^1-2tt 2 M 6 V- (4.4) 

To maximize the possible effect, we allow for the small deviations from the radial direction due to the finite size of 
the black hole. From simple geometry (see Fig. [2]), the particle will be moving at an average angle 

e bh (r) « — (4.5) 

r 

from the radial direction. Note that, for r >> rbh, almost all of the velocity is radial and Obh — > as r — > oo. The 
transverse component of the particle velocity in the black hole rest frame is 

v T ^ r -^, (4.6) 
r 



giving a radial velocity component of 



v R a 1 - 27 r 2 M 6 V - 2 -^k. (4.7) 



The radial velocity is approximately equal to vqm-, the velocity of the center-of-mass (CM) of two typical particles 
moving at an angle 9{,h either side of the radial direction. The corresponding Lorentz factor for the CM of the particles 
is 

At large r, this reduces to jcm ~ ^T^h/fA, as expected. The special relativistic transformation of velocities is 



v'=^^- v> = V -^^ (49) 
x 1-W v l-Vv x ' 1 ' 

where V is the relative velocity of the reference frames, taken to be in the x direction. Thus the relative velocity of 
the two particles in their CM frame is 

< e(K (l + ^V)" VS . (4-10) 

where the prime indicates the CM frame. This is close to 1 for radii less than /i -1 but very small at distances much 
larger than this. 

The total number density of emitted particles at a radius r in the black hole rest frame is given by Eq. (|3.ip . In 
the CM frame of the two particles, the number density at r is therefore 

»' W " W m 2 - !0" 3 ^ 3 (1 + TrVr 2 ) 172 ■ (4.11) 
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Hence the number of interactions per particle, integrated over all radii, is 

7V« / n'av' rel dt' = ( ^^dr. (4.12) 

J J 1CM 

Taking the cross-section in the CM frame to be the electron bremsstrahlung cross-section (the pair-production cross- 
section is similar), and allowing for the plasma mass correction (see Section III), we find 

10a 3 , A „, 



which gives 



10-V" 3 (1 + TT 2 m 2 e r 2 ) 1/2 a 3 (m 2 e + * M^r^dr 



for Tbh >> m e . The dominant contribution to the integral comes from the lower limit, here taken to be the radius of 
the black hole. Near this limit, the second term in the first brackets and the first term in the second brackets can be 
neglected, so we obtain 

POO 

W « / Wa 2 M bh r- 2 dr « 10a 2 << 1. (4.15) 



Hence emitted particles traveling in a similar radial direction will not interact with each other via bremsstrahlung or 
pair production, or scatter or lose energy to any significant degree, as they propagate away from the black hole. This 
is a consequence of the Lorentz factor of the emitted particles being so large in the black hole frame. 



C. Center-of-Momentum Frame of Interactions 



We now examine the Heckler assumption that the CM frame for most interactions is the rest frame of the black 
hole, so that (Jbrem is given by Eq. (|3.2|) in the black hole frame. As we have seen, this is not true for particles moving 
in similar radial directions. Here we investigate particles moving in different radial directions. 

Consider ultra-relativistic particles 1 and 2 of the same energy and rest mass, created near the black hole at the 
same time and subsequently moving radially outwards with speed v\ along rays separated by an angle 9. The Lorentz 
factor 7cm of the CM of the two particles relative to the rest frame of the black hole is 

7cm = [l-vl cos 2 - « sin- . (4.16) 



2 J V 2 , 

This is approximately independent of v\ (and hence Tbh) when 9 > 2/71. Equation (|4.16[) implies that 7cm — 100 
when 9 w 1°; 7c m = 10 when 9 w 11.5°; 7c m ~ 3 when 9 = 39°; and 7cM = 1 when 9 = 180°. Thus the CM 
frames for the interactions of particle 1 with the majority of particles streaming out from the black hole are not 
significantly Lorentz-boosted relative to the black hole. The region for which the interaction frames arc significantly 
Lorentz-boosted relative to the black hole is a cone of angle << ir/2 around the direction of particle 1. 

In general, the energies of particles 1 and 2 are not equal but have the Hawking distribution, given by Eq. (|2.1[) . 
and the corresponding formula for jcm is more complicated. Let us denote the ultra-relativistic speeds of particles 1 
and 2 by v\ = 1 — £\ and 1)2 = 1 — £2, respectively, where £1, £2 << 1, and let 9c mi be the angle between V\ and the 
velocity of the CM of particles 1 and 2, vqm- The special relativistic equations for relative velocity then give 

VCM « 1 ~( £l+£2 ) (4 17) 

((1 - cos 9) (^mi) + cos0cmi) - (£1 - e 2 cos9) ( s -^^) - e 2 cos9 CM1 



and 

sin6> CM i / £1 ((£1 + £2) - 2 v /£ie 2 cos6 



1/2 



sin6» I ( £l - e 2 f + 4£i£ 2 sin 2 £ 



(4.18) 



for 0° < 9 < 180°. Again 7cm = (l — v cm) IS approximately independent of V\ and v 2 (and thus Tbh) unless 9 
is extremely close to zero. If we assume E\ < E 2 , without loss of generality, and consider E 2 < 2E\, corresponding 
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to £2 < £i/4 (this more than adequately covers the majority of the Hawking distribution around E pea k), we have 
7cm < 3 when 9 > 42° and y C M < 10 when 9 > 12°. Similarly if E 2 < 4E X , then y C M < 3 when 9 > 50° and 
7cm < 10 when 9 > 15°. 

Hence we conclude that taking the CM frame of each interaction to be the black hole rest frame, as Heckler does, 
is a valid approximation for the interactions of any given particle with most other particles. To within a factor of 
O(l), using Ubrem in the black hole frame to describe the interactions of the particle is correct, except for an exclusion 
cone of angle 9 ex « O (0.1 — 1) radians around the direction of the particle's velocity for which the interactions are 
negligible because 7cm (9) >> 1. However, as we shall see below, this exclusion cone has important consequences 
which hinder the development of the photosphere. 

D. Transverse Distance 

Once the particle has traveled a distance d from the black hole, the transverse distance between the particle and 
the edge of its exclusion cone is 

xt « d6 ex « O (d) . (4.19) 

Because xt is a transverse distance, xt is the same in the particle and black hole frames. Thus once the particle has 
traveled a distance r\> rem f=a l/m e from the black hole, the distance to the nearest particle with which it could then 
interact is of order l/m e , provided the emitted particles are still traveling mainly radially. Because Ubrem dictates 
that two particles must be within ~ l/m e of each other to interact via brcmsstrahlung, few interactions can occur 
after the particles have propagated this distance. The relevant interactions for initiating photosphere formation can 
only be those which occur within rbrem of the black hole. This point is not included in the original Heckler model. 

As we shall demonstrate in Section IV. G, the electron coming out of a bremsstrahlung interaction is still off-shell at 
rbrem- Also Ubrem is truncated by the causality constraint and may be damped by off-shell (LPM-type) interactions, 
as we discuss in Sections IV. F and IV. H. These effects significantly decrease the capacity for photosphere formation. 

E. Deviations from the Radial Direction 

We have assumed in the exclusion cone argument that the emitted particles move out radially from the black hole. 
We now justify this. Were they to random walk out, an individual particle would have to undergo 




(4.20) 



scatterings to deviate by an angle 9 ex ss 0(0.1 — 1) from the radial direction, where <j) av is given by Eq. (|4.ip . If 
Tbh ~ 1 — 10 GeV, Nex ^ 10 7 scatterings are required for the particle to scatter outside the exclusion cone. Thus 
even if a photosphere can develop, particles will deviate little from the radial direction before and after photosphere 
formation. (As discussed in Section IV. L, this is incorrectly treated in the numerical work of Cline, Mostoslavsky, 
and Servant [l6|.) Therefore the constraint of Section IV. D set by the transverse distance to the edge of the exclusion 
cone will also apply after the onset of scattering. This confirms that the only relevant interactions can be those which 
occur within rbrem of the hole. 

Additionally if we consider the possible deviation from the radial direction due to the finite size of the black hole, 
Eq. (|4.5p implies that 9bh ~ 9 ex only when r ~ Tbh- This is approximately the radius at which the Hawking emitted 
particle first appears and, for Tbh » rn e , represents a tiny portion of the region within rbrem of the black hole. Hence 
the effect of black hole size on the propagation direction should also be negligible. 

F. Geometrical and Causality Considerations 

We now calculate the number of interactions Af each particle undergoes. The definition of the mean-free-path used 
in the Heckler model and given by Eq. (|3. 5|) is clearly wrong for geometrical reasons. Equation (|3 . 5|) applies when 
a particle propagates linearly through a target of uniform number density n in the configuration shown in Fig. [3] 
[381 ]. This is not geometrically analogous to the black hole situation for three reasons. First, the particles are moving 
radially outward in a spherical distribution from the black hole. Second, the number density n (r) varies with the 
distance of the moving particles from the black hole. That is, if we consider particle 1 interacting with an element 
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FIG. 3: The distribution of background particles for which the formula A = (nav)~ l describes the mean-free-path of an incident 
particle of speed v. The number density n is assumed to be uniform. In the actual black hole situation the density is not 
uniform and all particles are moving radially outwards in the black hole frame. 



of the flux around particle 2, n (r 2 ) does not remain constant as particles 1 and 2 travel outwards. Third, implicit in 
Eq. (|3.5[) is the assumption that the interacting particles have an infinite past and future, which is not the case for 
black hole emission. Consider particle 1 at radial position r x and particle 2 at r 2 with \r 2 \ < \ri\. Particle 2 was 
emitted by the black hole at a time At = |t , 2 | — | later than particle 1. The Hawking derivation of Eq. (|2.ip utilizes 
the solutions at infinity but does not tell us how a particle is evaporated: either the particle did not exist or it was 
hidden by the event horizon before emission. In either interpretation, the particle prior to its emission cannot interact 
with particles outside the event horizon. The interaction lifetime of the particle is therefore truncated. Because J\f in 
Eqs. (|3.5p and (|3.9|) is dominated by the interactions of particle 1 with the short-lived, high density particles close to 
the black hole, this truncated lifetime should have a significant effect on the Heckler scenario. 

An important consequence of the truncated lifetime is its implication for causality. From the classical point of 
view, i.e. ignoring temporarily the quantum smearing of particle space-time position and energy-momentum due to 
the Heiscnberg Uncertainty Principle and the quantum nature of interactions, two particles can interact only if they 
are causally connected or, cquivalcntly, if a light signal can be transmitted between them. Classically each particle is 
point-like and emitted by the black hole at a precise time. Without loss of generality, let us assume particle 1, emitted 
at time t\ in the black hole rest frame, is the scattering particle in the bremsstrahlung interaction and particle 2, 
emitted at t 2 = t\ + At, is the scattered particle. The time difference At may be positive or negative. The most 
appropriate frame for analyzing the causality constraint is the frame of the scattering particle because the Coulomb 
field of the scattering particle has an undistorted e/r 2 distribution in its own frame. Let us consider a light signal 
which travels a distance D from the scattering particle in the scattering particle's frame and reaches the scattered 
particle. The form of the invariant <7b rem , Eq. (|3.2|) . implies that, if the particle histories were not truncated, the 
probability for the scattered particle to undergo a bremsstrahlung event would be ~ a 3 if D < l/m e , but negligible 
if D >> l/m e . (The truncation of the particle histories will further suppress the probability within D < l/m e j39j.) 
If At > and particle 1 emits a light signal which reaches particle 2 just as particle 2 is emitted by the black hole, 
one can readily calculate [39[ that 

D = vV - 1 At = jvAt (4.21) 

where v is the velocity of the particles in the black hole frame. If, alternatively, At < and particle 1, as soon as it 
is evaporated by the black hole, emits a light signal which then reaches particle 2, a slightly longer calculation (39j 
shows that 



D = 7'i> 



1 — v cos 9 
1-v 



|Ai|>7«|Ai| (4.22) 



where 9 is the angle between the particles in the black hole frame. Thus in cither case, for the D < l/m e constraint 
to be met, the difference in the emission times must be 

At < At c = — !— a — (4.23) 

jvm e 7m e 
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for relativistic particles. However, from Section II, the average time between successive emissions of electrons and 
positrons directly evaporated by the black hole is 

*.-#5-. (4.24) 

^peak 

Using 7 f=a E pea k/m e , we see that 

At c « At e (4.25) 

for all Tbh >> fn e . Thus the classical causality condition (|4.23[) is only satisfied by a negligible fraction of the emitted 
particles when Tbh » m e . The same conclusion applies for QCD particles when T\,n » A-qcd- Because the emission 
of a given species is damped near its rest mass threshold |28| , so increasing At e , the classical causality condition is 
unlikely to be satisfied even near particle rest mass thresholds. (Also recall that Obrem decreases faster with decreasing 
E in the non-relativistic limit.) At any Tbh, if a- photosphere is to form, we additionally require the causality constraint 
to be met within a sphere of radius Tbrem ~ l/ TO e around the black hole. This only further tightens the At constraint. 

Quantum mechanical effects relax the above classical constraint only slightly. By the Heisenberg Uncertainty 
Principle, the emission of each electron from the black hole is smeared out over a time At m {pfm e ) « At e , 
leading to it having a spatial spread in the direction of emission of Ar as (^ym e )~ . Thus when the electron is a 
distance r from the black hole, we can visualize it quantum mechanically as an eggshell segment of thickness ('jm e ) , 
radial curvature r and tangential width growing to l/m e , whose probability of overlapping another emitted electron 
is essentially given by the classical constraint above. Additionally, the virtual photon exchanged between the two 
bremsstrahlung- interacting electrons, which typically transfers 3- momentum |fe 7 »| ~ m e and negligible energy and 
has a nonzero mass (miL as — m^), cannot be localized within a length less than l/|fe 7 »| sa l/m e . It can extend over 
a spacelike interval but spacelike intervals in the propagator are strongly suppressed. However, the wave speed at 
which the field of one electron propagates to the other electron in the Dirac equation remains c. Therefore the above 
classical causality constraint essentially always applies to the quantum interaction. 

Equation (|3.4[) then must be replaced with an expression which more accurately describes the geometry of the 
interaction and accounts for the shortened histories of the interacting particles. As we discuss in Ref. [39(, the causality 
constraint and truncated histories also decrease the momentum transferred in each interaction, decreasing the effective 
&brem WC H below the already small (Tbrem oc a 3 cross-section of Eq. (|3.2[) . Even at high Tbh, the average number of 
interactions per particle is bounded by At c / At e as 1CU 2 and cannot reach Af = 1. If an emitted particle does scatter 
once, we now show that the distance to complete that first scattering is such that subsequent complete scatterings 
by the particle within Tbrem of the black hole are not possible. Thus it is improbable for particles to experience more 
than one bremsstrahlung interaction. This by itself is sufficient to prevent QED photosphere formation. 



G. Scale for Bremsstrahlung Interaction 

In the original Heckler model the distance needed for the completion of a bremsstrahlung interaction, i.e. the 
distance required for the formation of the final on-shcll electron, is assumed to be df orm as l/m e in the CM frame of 
the interaction. This is not correct and must be replaced by the much larger distance df orrn as E/m>g, as we can see 
in two ways. 

First, the 4- momentum of the off-shell electron is — m \ + 2pi-k where p x and k are the 4-momenta of the final 
on-shell electron and photon, respectively (see Fig. 1). From Eqs. (|4.1[) and (|4.2p . the average energies of the final 
electron and photon are ~ E/2 and the average angle between them is cf) av as m e /2E, so we have p^ — ml as ml/2. 
Thus the average lifetime of the off-shell electron in its own frame is ~ l/m e , which corresponds to ~ E/m 2 e in the 
CM frame of the interaction. 

Second, following Ref. consider the final electron and photon as they are being created. Because they are 
created with a small separation angle, their wave packets will initially overlap for some distance as they propagate. 
This distance is the formation distance. The transverse momentum components of the final electron and photon 
relative to the momentum of the initial electron p 1 are both px ~ 4>avE. By the Heisenberg Uncertainty Principle 
these transverse momenta can be resolved only when the transverse spatial separation is at least 1/pt- Since the 
opening angle between them is <f> a v, this occurs after they have traveled a radial distance (j>T<pav) ~ Ejm\ in the 
CM frame of the interaction. Thus again the formation distance must be df orm as E/m\. 

As we have already noted, the consequence of this is that the only relevant interactions for photosphere creation 
can be those which occur within Tbrem ~ l/m e of the black hole. Certainly df orm >> l/m e for Tbh >> m e and 
so, on completion of the first bremsstrahlung interaction, the products can no longer significantly interact with other 
particles. In the extremely rare event that an electron underg oes more than one scattering within T^ rem of the black 
hole, the subsequent scattering occurs when the electron is off-shell, not on-shell. 
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H. Off-Shell Interactions and the LPM effect 



In the general setting, if an off-shell particle undergoes multiple scattering over the distance df orrn (the formation 
length of the photon and final electron in the bremsstrahlung interaction), Landau, Pomeranchuk and Migdal have 
shown that the bremsstrahlung cross-section is decreased by a numerical factor of up to 0(1) (see [4l| ; for a recent 
review see Ref. 1431 ). This is known as the LPM effect and has been confirmed experimentally, at least in the soft 
photon regime [431 ] . (The effect is greatest in this regime because df orm increases as to decreases |44|.) The suppression 
of the bremsstrahlung cross-section is due to multiple scattering breaking the coherence of the final state as the photon 
is forming [45[. The suppression increases with electron energy and the number of off-shell scatterings or, equivalently, 
the number density. The LPM effect modifies both QED and QCD bremsstrahlung. Recent theoretical treatments 
include elastic (Coulomb) and inelastic (bremsstrahlung) scatterings of the off-shell electron [Z(| EH . 

Applying the LPM effect to the black hole scenario, we would expect multiple scatterings of the off-shell electron to 
decrease the effective bremsstrahlung cross-section. Since the off-shell electron requires a distance df orm to decay into 
the final on-shcll electron and photon, even if J\f > 1 within rbrem l/m e , each particle generally undergoes at most 
one complete on-shell scattering as it streams away from the hole. Recall that the probability of undergoing even one 
complete scattering is significantly damped by the causality constraint at >> m e . Incorporating the precise effect 
in the black hole scenario requires calculating the LPM effect from first principles because the black hole scenario differs 
from other LPM scenarios in at least three ways: (i) the black hole off-shell interactions are bremsstrahlung (not just 
Coulomb scatterings) with average momentum transfer similar to that of the initial bremsstrahlung interaction; (ii) 
the off-shell interactions and the initial interaction occur essentially simultaneously, not sequentially, and their range 
is greater than the distance between each off-shell interaction; and (iii) the non-uniform number density, truncated 
particle histories and energy distributions of the emitted particles are specific to the black hole case. 

Proceeding qualitatively, we expect the following outcomes for black hole QED interactions. For bremsstrahlung, we 
have r^ rem sa l/m e << df orm , so multiple scatterings, if they occur, can perturb only a small segment of the complete 
interaction of a typical particle. Because df orm » r\, remi any scatterings affect the particle essentially simultaneously 
within a distance r < rtrem of the black hole. (The actual region will be smaller than r^ rem because of the causality 
constraint.) Therefore the particle can be thought of as interacting with a region within r < rbrem as a whole. This 
corresponds to the single radiator or thin target LPM limit [42|, Es l48l . The probability for interaction in this 
limit is given by the total scattering angle within the r < rbrem region [42|. If we approximate this angle by that 
appropriate for a random walk of successive scatterings, 9ms ~ VNs<fiav where Ms is the number of scatterings the 
particle undergoes within r < rbrem of the black hole, then 9ms > 4>av for A/jg > 1. From the numerical simulations 
of Haug 3J|, dtJbrem / d9 falls off roughly exponentially at angles above <p av and the momentum transfer is increased, 
although the energies of the final photon and electron stay of order E. Hence dubrem {9ms) /dO < dobrem {4>av) /dO 
and the modified cross-section is always less than that for the interaction without off-shell scatterin g. A dditionally, if 
the bremsstrahlung photon pair-creates before it is completely formed, <Jbrem is further suppressed [421 ] . 

Thus the net effect of multiple off-shell QED scatterings is to decrease the probability of any particle producing a 
bremsstrahlung photon. Even if J\f > 1, each particle requires a distance df orm » rbrem to complete the interaction 
and so is unlikely to undergo more than one completed bremsstrahlung interaction as it streams away from the black 
hole. Because each particle rarely undergoes more than one completed interaction, the number density around the 
black hole can increase by only a small fraction (at most 50%) even if Af > 1. Hence we conclude that a QED 
photosphere will not develop. 



I. The Plasma Mass Correction 



In the original Heckler model, it was argued that on scales less than rb re m ~ 1 /m e the QED bremsstrahlung cross- 
section should be suppressed by Coulomb scattering off other electrons around the black hole. (However, Heckler 
did not include the off-shell bremsstrahlung interactions discussed above.) To account for the Coulomb scattering, 
Refs. [l3|, [lj| took the simplified approach of defining the effective electron mass m e to include the fermion self-energy 
for a finite temperature bath. In the expression for abrem, rn e was augmented in the interaction rest frame by the 
thermal plasma mass given by Eq. (|3.3|) . which describes the screening of the electron by other emitted electrons. As 
Heckler pointed out, the Hawking emission is not a true plasma. Because the emitted particles flow isotropically out 
of the point-like black hole, the particle distribution cannot be transformed into a frame in which the particles exhibit 
a Maxwell-Boltzmann momentum distribution. 

As we have seen, the causality constraint prevents emitted electrons interacting significantly within rbrem ~ 1/ m e 
of the black hole. This should apply to Coulomb scattering as well as bremsstrahlung interactions. Also, since 
each electron requires a distance df orm >> rbrem to complete a bremsstrahlung interaction, n(r) cannot increase 
substantially within rbrem of the black hole due to interactions. Thus plasma mass considerations should be irrelevant 
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within r brem w l/m e for T 6ft » m e . 

J. Other Bremsstrahlung Processes 

Although the Heckler photosphere scenario relies on the O (a 3 ) two-body bremsstrahlung interaction between two 
charged particles emitted by the black hole, the emitted particles may experience other QED bremsstrahlung processes 
which we investigate in detail in the accompanying paper (39l |. First, there should be the O (a) inner bremsstrahlung 
generated by the charged particle as it escapes the black hole and, from the point of view of a distant observer, has 
its velocity changed from being essentially zero at the black hole to its asymptotic value. This inner bremsstrahlung 
emission will contribute about 10% of the total photon power from an Mbh = 5 x 10 14 g black hole and has an 
approximately energy-independent power spectrum which dominates over the directly emitted Hawking photons at 
low energies J39|. Secondly, there should be the O (a 2 ) bremsstrahlung emission from the scattering of the charged 
emitted particles by the stochastic electromagnetic field of the black hole, and another O (a 3 ) contribution from the 
scattering of the inner bremsstrahlung photon off an ambient charged particle. In fact, since the inner bremsstrahlung 
dominates all bremsstrahlung processes, the other bremsstrahlung processes can be simply regarded as small O (a) 
and O (a 2 ) corrections to the inner bremsstrahlung. The total bremsstrahlung emission is then essentially determined 
by the asymptotic momentum distribution of the charged particles, which to lowest order in a is determined by the 
Hawking emission formula Eq. (|2.1| and to higher orders is modified slightly by the interactions between the emitted 
charged particles. In all cases, however, the outgoing photon and electron from any bremsstrahlung interaction require 
a formation distance df orm ~ j/m e . Thus none of these bremsstrahlung interactions can produce a QED photosphere 
around a Tbh >> n^e black hole. 

K. Extension to Pair-Production and QCD Effects 

So far, we have discussed QED bremsstrahlung interactions but the above considerations also apply to QED pair 
production e + 7^e + e + + e~. The pair-production cross-section is obtained from Obrem by rotating the Feynman 
diagrams to depict an incident photon rather than an outgoing photon and using the substitution rules (34[. In the 
relativistic limit, cr pa , r has the same form as Obrem [35( ] . The distribution of the outgoing e + and e~ is strongly peaked 
in the direction of the initial photon and each carries on average approximately half the initial photon energy [3^ |. 
The formation length of the e + e~ pair is also df orm ~ w/2m 2 >> l/m e when Tbh >> l/m e [441 ] . As before, any 
multiple off-shell scattering will suppress a va ir within the distance df orm when Tbh » fn e . Thus each particle can 
undergo at most one complete pair-production event as it streams away from the black hole and the probability of 
that interaction is strongly suppressed by the causality constraint. Again no QED photosphere can develop. 

The above discussion also applies to QCD gluon-bremsstrahlung and QCD pair-production interactions of the 
emitted particles. (Note here we are referring to the QCD interactions between particles in different QCD jets, 
not between particles in the same QCD jet. The interactions between particles in the same jet are included in 
the fragmentation and hadronization modeling of Ref. [12].) To derive the corresponding QCD quantities, m e is 
replaced by the quark mass m q or Aqcd, and a by the strong coupling constant a s . This leads to Eq. (|3.14p in the 
relativistic case, neglecting color factors of O (1). The average time between QCD particle emissions by the black 
hole is AtgcD ~ 20/Ep ea k when Tbh >> Aqcd ■ The QCD formation scale d^J^ « E/Aq CD is also the scale at 
which the Hawking emission hadronizes. Applying the same arguments as above, when Tbh » Aqcd, the causality 
constraint implies that the emitted QCD particles have little chance of interacting within r2. e ™ ~ 1 /Aqcd of the black 
hole and any particles that do interact can undergo at most one complete QCD bremsstrahlung interaction within 
^form >;> IMqcd- Thus no QCD photosphere (also known as a chromosphere) should form when Tbh >> Aqcd- 
Because of the time separation between subsequent QCD jets and the LPM effect, nor should a photosphere form 
from interaction of the hadronization products which appear at d^J^. 

QCD phenomenology and experiments also strongly argue against a QCD photosphere developing when Tbh (or 
more specifically the peak energy of the directly emitted flux) passes through values that are roughly Aqcd- Because 

the Heckler estimate, T2-f t D ~ Aqcd / &V 2 , is close enough to Aqcd that high energy formulae may be inappropriate, 
the question arises of whether a QCD photosphere could form from non-relativistic effects at these temperatures. If 
the emission were non-relativistic, the causality constraint would be weakened. Also at low energies the effective a s 
is observed to increase from its high energy limit, a s = \2 ir/ { (33 — 2n/) In (E 2 /Aq CD ) } where n/ is the number of 
relevant quark species, and may approach 7r as E — > Aqcd [4^[50l.[5ll. |52|. On the other hand, it is not clear whether, 
even if multiple scatterings were to occur around Tbh ~ Aqcd, they would affect the observable spectra, because 
the Hawking emission is strongly damped and the multiplicity per QCD jet is very low around the Aqcd threshold 
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[12|. However, we deduce from both the QCD phenomenology and experimental perspectives that when the black 
hole begins to emit asymptotically free quarks, these quarks are relativistic. Therefore no QCD photosphere should 
develop around Tbh ~ ^qcd- 

The emission picture which is consistent with QCD phenomenology and threshold accelerator experiments is as 
follows. When Tbh (or more precisely E pea k as 2.81Tj,/j for s = 0) is about the pion rest mass (m n ± as 140 MeV 
and m^o as 135 MeV), the black hole should start emitting pions directly, with the grey-body factors and degrees 
of freedom appropriate for massive s = particles modified by the finite-size structure of the pion. Once the scale 
of the black hole emitting region becomes small compared with the length below which the constituent (valence) u 
and d quarks in the pions become asymptotically free, the black hole should begin emitting the quarks individually. 
These quarks then hadronize into pions on the proper distance scale 1 / Aqc d m the frame of the emitted quark. The 
energy at which this happens is sufficiently above the pion rest mass that the pion is relativistic. Hence the pion's 
constituent quarks should also be regarded as relativistic. Furthermore, during the asymptotic freedom regime before 
the quark hadronizes, the quark masses which are used in the QCD Lagrangian are the much lighter quark current 
masses, not the quark constituent masses. (The quark current masses are m u as (1.5 — 3) MeV, as (3 — 7) MeV 
and m s as (70 — 120) MeV, whereas the quark constituent masses are m Ut d as 0.31 GeV and m s as 0.48 GeV. For the 
heavy c, b and t quarks, the current masses are approximately equal to the constituent masses and lie above 1 GeV 
[53l ]. The quark constituent masses are the relevant masses as the quarks hadronize and afterward.) Thus the quarks 
are ultra-relativistic during their asymptotic freedom phase and then gain larger effective masses, but still remain 
relativistic, when they hadronize. 

The hadronization process itself at these energies also cannot produce a significant increase in the number of final 
particle states. What the quarks hadronize into is dictated by energy and quantum conservation laws. Hence, around 
~ h-QCD, the quarks can hadronize only into pions (or the few light meson states which decay into pions or 
mimic ir° decay into photons) and the m ultip licity of the final pions per initial Hawking emitted quark must be very 
low (one or two pions per initial quark) [531 ] - Emission of e + e~ or photons by these quarks or pions is covered by 
our earlier analysis of QED interactions and so is negligible, as are any weak interactions. Gluon bremsstrahlung 
by quarks at these energies is also inconsequential. The gluon has to obey the QCD confinement (no-free-color) and 
quantum conservation rules. Thus a gluon must have an energy of about Aqcd or greater to produce free final state 
particles, the lowest energy option being a pion pair. Because the total momentum in the frame defined by the black 
hole and the Hawking emitted quark must be conserved and there are no available very light final states for the gluon 
or the quark, an asymptotically free Hawking emitted quark of energy E ~ Aqcd cannot be slowed down by gluon 
bremsstrahlung, except by the gluons involved when that quark eventually hadronizes. 

The above description is completely consistent with accelerator experiments. The cross-section for e + + e~ — > pions 
has been measured at center-of-mass energies from close to the 2m 1T threshold to 1 GeV and remains smooth over that 
region (see Ref. [54[ for a recent review). No 'dip' or other structure in the data is observed around center-of-mass 
energies of Aqc d where the behavior changes from direct production of pions to the production of final state pions 
via an intermediate asymptotically free quark state. 

To further visualize the black hole quark emission and hadronization process around Tbh ~ ^-qcd and at higher 
temperatures, one can use the chromoelectric flux-tube model (also known as the QCD string model) that successfully 
describes hadronization in accelerators (see [55[, (56|; for recent reviews see Refs. [UHll). In this interpretation a 
quark is emitted from the black hole dragging behind it a tube of color field whose other end is attached to the black 
hole. The string tension cta of the flux-tube is constant, so that its potential energy a^x grows linearly with the length 
of the tube x. When the flux-tube has been stretched to the length x as y / 2n/a\ as 1/Aqcd at which the potential 
energy matches the energy required to create a new qq pair out of the vacuum (including the transverse momentum of 
the qq pair implied by the Hciscnbcrg Uncertainty Principle applied to the flux-tube width), the tube breaks in two, 
producing a q and a q attached to the two new ends where the tube broke (57| . To conserve energy, the stretching and 
qq creation is accompanied by a corresponding deceleration of the lead quark. This stretching and breaking sequence 
continues until the energy of the lead quark has been expended producing new qq pairs. The quarks and antiquarks 
created from the flux-tube match up to hadronize into bound color-neutral states. The whole process is completed 
over the distance scale E/Aq CD in the black hole frame, where E is the energy of the initial quark. The flux-tube 
model naturally gives a multiplicity growth of final states per initial quark that is proportional to \aE, as observed 
in accelerators. 

That the asymptotically free quarks arc relativistic, rather than non-rclativistic, when first emitted by the black hole 
around Tbh ~ Aqcd can additionally be seen from the thermodynamical 'detailed balance' derivation of the Hawking 
flux. When a black hole is in quasithcrmal equilibrium at a given instant with a surrounding thermal bath of the same 
temperature, the absorption rate of a given species equals its emission rate. If a pion incident upon a black hole has 
total energy insufficient for the pion to break up into subparticles with one subparticle going down the hole and one 
escaping, detailed balance gives the Hawking emission rate of direct outgoing pions from the capture cross-section of 
the incoming pions. When the pion size is comparable to the cross-section for a point-like scalar particle of the same 
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mass and energy as the pion, the pion capture cross-section should be similar to that of the corresponding point-like 
scalar particle. (This cross-section should also apply when the black hole itself is much smaller than the pion if the 
pion is moving sufficiently slowly that the capture cross-section is much larger than the hole itself. In this case, once 
the pion falls to a distance from the hole comparable to its size, it will be nearly inevitable that part of the pion 
will fall into the hole, with the rest having insufficient energy to escape even though the pion may initially appear to 
have too large an intrinsic size to be captured by the smaller hole.) If the black hole size is significantly smaller than 
the size of the pion and the incoming pion is moving rclativistically (that is, its kinetic energy is significantly greater 
than its rest mass), the black hole can capture an individual constituent quark of the pion, with the other constituent 
escaping capture or being captured after a significant time interval. This inelastic scattering of the pion by the black 
hole could also lead to the production of other particles, depending on the available energy. One could visualize the 
capture process as the glue string between the captured quark and the constituent antiquark of the original pion 
breaking in two, with a new antiquark forming on the end of the string attached to the quark that fell into the hole 
and a new quark forming at the end of the string attached to the uncaptured antiquark. The new antiquark attached 
to the glue string going into the hole could be captured by the hole, while the new quark attached to the external glue 
string could escape along with the original antiquark. For incident pions of high enough energy, more complicated 
inelastic processes could lead to greater multiplicity and number of species in the escaping and captured states. The 
reverse process of the absorption of constituent quarks corresponds to the Hawking emission of individual relativistic 
quarks that hadronize into one or more pions and/or other mesons and baryons. 

From our above discussion it is clear that the QCD production around a T^n ~ Aqcd black hole is not analogous 
to the creation ofquark-gluon plasma at the Relativistic Heavy Ion Collider (RHIC) (for a recent review of RHIC 
results, see Ref. |50|). At RHIC the collision of heavy 197 Au + 197 Au nuclei at center-of-mass energies of 130 GeV 
and 200 GeV per nucleon generates a strongly coupled, high entropy quark-gluon plasma which rapidly achieves a 
thermal spectrum and then expands iscntropically until the particles hadronize and decouple by a time 10 — 20 fm/c. 
Presently it is not known whether the plasma achieves its thermal spectrum of states by interactions or is initially 
created with that distribution of states jsjj. At the moment of collision, the overlapping 197 Au nuclei arc highly 
Lorentz-contracted in the center-of-momentum frame, implying an overlap time of 0.2 fm/c and an initial thermalized 
energy density of about 15 — 3000 GeV fnV 3 . The initial energy density estimate varies with the interpretative model 
but in all cases is far greater than the 1 GeV fm -3 threshold for quark-gluon plasma formation or the density within 
a nucleon 0.5 GeV fm -3 , and remains so for up to 5 fm/c. The number of charged particles arriving at the detectors 
(about 5000) indicates that approximately 20 high energy particles are created per initial nucleon but the dependence 
of the charged particle multiplicity on the available energy per colliding nucleon still matches that seen in e + + e~ and 
p+p accelerator collisions [601 ] . The hadronic multiplicities are determined by the nonzero baryonic chemical potential 
(baryon/anti-baryon asymmetry) of the initial nuclei. Notably, a reduction of a factor of about 5 in the number of high 
transverse momentum hadrons is observed. This phenomenon is known as jet-quenching, although the high transverse 
momentum jets are not truly quenched but spread over a wider angle with lower average energy per particle. The 
hadronic energy loss implies significant QCD interaction with the ambient plasma, but photons are observed to travel 
through the plasma unimpeded. The jet-quenching can be modeled by assuming that the collision initially creates a 
gluon number distribution of dN g /dy ~ 1000, where y is the rapidity defined by E = yJ m? + pip cosh y for a particle 
of energy E — jm and transverse momentum px- (Integrating over 1 < coshy < 7, this gives the number of gluons 
per nucleon to be of order 10 — 10 2 ). Application of the Heiscnbcrg Uncertainty Principle indicates that the nuclei 
are gluon-saturated on the scale of the collision even before impact: that is, at 200 GeV per nucleon the two gold 
nuclei see each other as dense objects, not ensembles of distinct partons, because the atomic number (parton density) 
and transverse momenta are such that the partons overlap considerably in the transverse plane [6l| . All of the above 
effects attributable to a strongly-coupled quark-gluon plasma are experimentally confirmed to be suppressed at lower 
collision energies and in deuteron-gold collisions [62j ■ 

In the Tbh ~ ^qcd black hole case, however, in contrast to RHIC, the conditions for the simultaneous production 
of an ultra- high density of QCD particles are not met. A T^h ~ A-qcd black hole emits particles by the Hawking 
process sequentially and relativistically with significant time separation between each particle emission. Thus, as we 
have seen in Section IV. F, negligibly few of the Hawking emitted quarks are in sufficient causal contact to interact 
with each other within 1/Aqcd of the black hole. Furthermore, as discussed above, conservation of energy severely 
limits the number of final states that can be created per initial quark around T^n ~ ^qcd- In the black hole case, only 
an energy E ~ Aqcd is available per Hawking emitted quark, in contrast to the 200 GeV per nucleon available at 
RHIC. (Additional differences include the facts that the corresponding baryonic chemical potential is fis = 0, because 
the black hole emits equal numbers of quarks and antiquarks, and that the black hole particles arc not initially gluon- 
saturated.) Hence, as given by Eqs. (|2.1ip - (|2 . 15[) . the greatest achievable density around a Tbh ~ Aqcd black hole is 
well below 1 GeV fm -3 . The closest experimental analogy to the T^h ~ A-qcd black hole situation is the generation of 
pions and QCD jets in e + + e~ accelerator collisions and in hadronic collisions at much lower center-of-mass energies 
per nucleon. 
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In summary we conclude that the time separation between emissions and the limited average energy available per 
Hawking emitted particle prevent the formation of a QCD photosphere around a Tbh ~ Aqqd black hole, and the 

causality and dy^^ >> 1/A.qcd constraints prevent QCD photosphere formation at Tbh » Aqcd- We also note 
that the above description of the QCD threshold emission processes should lead to more precise numerical modeling 
of the astrophysical particle spectra produced by Hawking radiation around the Aqcd threshold than has previously 
been published. 

L. Comments on Cline, Mostoslavsky, and Servant Simulation of Heckler Photosphere 

The original Heckler scenario was numerically modelled by J. Cline, Mostoslavsky, and Servant [Hj], who simulated 
photosphere development usin g th e test particle method to solve the Boltzmann equation for propagation of plasma 
particles undergoing collisions [lfj]. They found that QED photosphere production set in at T cr u R* 50 GeV (compa- 
rable to Heckler's value of 45 GeV) but extended over a much smaller region: it started at n n ~ I0 4 rbh (compared 
to Heckler's n„ sa lO 9 ^) and ended at r ou t ~ l/m e where the trajectories of particles within the interaction range 
rapidly become parallel, as we saw in Section IV. D. For Tf,h > 50 TeV, their photosphere region was extended by 
Compton scattering effects. However, Cline, Mostoslavsky, and Servant use ac = 2Tra 2 \n(E/m e )/(m e E), whereas 
the correct relativistic Klein-Nishina cross-section [63[ for a photon of energy E Compton scattering off an electron 
of energy E in the black hole rest frame is ac ~ 2ira 2 \n(E /m e )/ 'E 2 . They also found that the fluid description is of 
questionable validity and that it led to an average photon energy at the outer photosphere edge which differed from 
Heckler's value, even decreasing as Tbh increases. The differences are more marked in the QCD case, partly because 
of the form of the QCD plasma mass correction used by Heckler. They found that a QCD photosphere forms for 
Tbh > T? r ? t D « 175 MeV and extends to r out « 1/Aq C d- 

The Cline, Mostoslavsky, and Servant work is elucidating but, because it takes the same approach as Ref. all of 
the points we have made above apply to it. Most importantly, the authors neglect the causality and df orm constraints 
which limit each particle to at most one on-shcll brcmsstrahlung interaction within r ou t f=a l/m e or r ou t w 1/Aqcd 
and make even that highly unlikely. An additional problem is that they generate the bremsstrahlung photons with 
random energy < oj < E and random angle < cf> < 2n. We have seen that the true photon angle is very small 
(j) av ?s m e /2E and deviates little from the direction of the initial electron. Randomizing the angle would exaggerate 
photosphere development, and this may explain why Cline, Mostoslavsky, and Servant find T cr i t to be lower than 
Heckler's original analytic estimate. Their numerical method should therefore be modified to include all of our previous 
points in Section IV and could be employed to investigate the consequences of the rare bremsstrahlung interactions. 

M. Comments on Kapusta and Daghigh Models 

Kapusta [13, EH and Daghigh and Kapusta H3, HH have analytically and numerically investigated the hydro- 
dynamic expansion of matter around a black hole in the original Heckler photosphere formation scenario. In their 
treatment, the Tbh > 100 GeV emission is assumed to be a quasi-stationary shell of plasma expanding relativistically 
away from the hole. The plasma is kept in local thermal equilibrium by the QED and QCD brcmsstrahlung and pair- 
production interactions of the original Heckler model. 'Local thermal equilibrium' in this context is taken to mean that 
the mean-free-path, or thermalizing scale, of the plasma particles is less than the scale over which the plasma temper- 
ature is changing They find that the quasi-stationary photosphere shell extends out to 10 -10 (Tb^/l TeV) 2 cm. 
The expected observable photon signal is then calculated from the luminosity of the photosphere surface. 

The results of Daghigh and Kapusta are not supported by our analysis. As we have seen in Section IV, the 
interactions within r < rbrem of the black hole cannot significantly increase the number density n (r) , as required for 
photosphere formation. Their approach also neglects the causality and df orm constraints within r < rbrem- Once a 
particle reaches rbrem-, we have shown that the geometry is such that it has little probability of interacting with other 
emitted particles. Thus, contrary to the claim of Daghigh and Kapusta, there should be no change in the particle 
energy after the particle reaches r w r^em and no persistent QED or QCD photosphere or locally thermalized plasma 
should develop. 

V. CONSEQUENCES FOR OBSERVABLE SPECTRA 

In Section IV we have shown that a Heckler QED photosphere will not form around a Hawking black hole. Similar 
arguments and additional QCD phcnomcnological considerations prevent a QCD photosphere developing. Interactions 
which occur within Tb rern of the black hole determine whether a Heckler photosphere can arise. These interactions 
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are extremely rare and any multiple scatterings involve off-shell particles. In contrast, the original Heckler model 
postulated that each particle on average undergoes successive on-shell bremsstrahlung interactions out to a much 
greater radius. This led to an increase in the expected flux at low energies and a steepening of the spectra at 
high energies. Including the off-shell scatterings, however, suppresses the effective bremsstrahlung cross-section and 
permits at most one completed bremsstrahlung interaction per particle. This one net interaction can only slightly 
modify the observable emission spectra, and not significantly as would a photosphere. Our remarks apply equally to 
pair production. Also of crucial importance to determining the full spectra is understanding the emission processes 
around T bh ~ A QC d- 

If each emitted particle produces at most one bremsstrahlung particle or undergoes pair-production once as it 
streams away from the black hole, there will be little effect on the detectability of an individual black hole or a 
background of black holes. The number of relativistically emitted particles of the relevant species in the relevant 
part of the spectra could be increased by at most a factor of 3/2 and the average energy decreased by at most 1/3. 
Furthermore, given the low likelihood of even one net interaction, the causality constraint discussed in Section IV. F 
implies that any change to the spectra should be several orders of magnitude weaker than this. As we show in 
Ref. [39l | , inclusion of the causality constraint in an impact-parameter-based analysis of the two-body bremsstrahlung 
interaction around a Tbh >> Tn e black hole leads to an upper limit of 10 -7 on the total fraction of energy emitted 
in two-body bremsstrahlung photons. Charged and/or colored particles (e^ : /i^, r^, q^, g) and photons which 
are directly emitted by the black hole could participate in the rare bremsstrahlung or pair-production interactions. 
Additionally, almost all of the 7, p, p, e + , e~ , v and V flux from a Tbh <^ 1 GeV black hole comes from the decays of 
directly emitted particles. Thus all observable spectra, except the highest energy vv flux which is directly emitted, 
may be very slightly modified by interactions. As we discuss in Ref. [39] . at low energies the photon spectrum from 
an individual black hole will be enhanced by two other O (a) and O (a 2 ) bremsstrahlung processes. In contrast, the 
low energy particles in the spectra from a PBH background were originally emitted with higher energies and then 
cosmologically redshifted to low energies [64| . and so the background spectra can be only slightly modified at low 
energies (by the rare interactions near the black hole). In summary, the observable spectra from an individual black 
hole or PBH background remain those of the standard Hawking picture. The time-integrated background spectra 
above E w 300 MeV should maintain their original dN/dE cx E~ 3 form [l(| with little amplitude change, in contrast 
to the E~ A slope predicted by the Heckler model [3]. The limits on Q p bh derived at E f« 100 MeV using the 7, p, p, 
e+, e~, v and V background spectra [l(| will not be significantly affected. 

We also note that the high energy limits on the detectability of individual holes, which are derived at E >> Aqqd, 
are not sensitive to behaviour around Tbh ~ ^-qcd but the limits on Cl p bh, which are derived at 100 MeV, may 
be. This is because at E m 100 MeV most particles in the integrated spectra from a PBH background are QCD 
jet products from Tbh <] ^qcd black holes emitting in the present era [64] , The Galactic antiproton limit is less 
influenced by the behaviour around Tbh ~ ^qcd because most of the postgalactic PBH antiprotons are generated in 
somewhat higher energy QCD jets, due to the greater rest mass of the antiproton [64[. In future work a more precise 
numerical modeling of the spectra at E « 100 MeV using our insights of Section IV. K should be pursued. 

VI. OTHER PHOTOSPHERE SCENARIOS 
A. Belyanin, Kocharovsky, and Kocharovsky Model 

Belyanin, Kocharovsky, and Kocharovsky [22j have proposed an unrelated photosphere model, in which charged 
particles from the black hole — predominantly pions from the decay of Hawking radiated quarks and gluons - 
form a well-defined 'collisionlcss' plasma at r » Tbh- The hydrodynamical behavior of this plasma is assumed to 
be maintained by collective wave-particle interactions and/or a self- induced turbulent magnetic field, both of which 
the authors postulate can arise from anisotropies in the charged particle distributions. Particle-particle collisions 
themselves are insufficient to maintain this hydrodynamical regime. Belyanin, Kocharovsky, and Kocharovsky claim 
that there is sufficient time for plasma wave turbulence to develop as the particles stream away from a high temperature 
black hole. Explicitly, they find that, when Tbh ^ 10 GeV, the size of the plasma (taken to be the distance from the 
black hole) is greater than the Debye radius ro (effectively the inverse of the Langmuir oscillation frequency — or 
plasma frequency generated by the charged particle separation) , and that ro is greater than the average separation 
between charged particles. 

The authors apply magneto-hydrodynamical equations for a plasma magnetic field in the case where the magnetic 
field is maximal, i.e. where the energy is equipartitioned between particle kinetic energy and the magnetic field. They 
argue that, for Tbh ^ 10 GeV, synchrotron radiation and electromagnetic cascading in the self-induced magnetic field 
will produce an observable MeV 7-ray burst lasting 10 _1 — 10 3 s. They envision three possible sources for the induced 
turbulent magnetic field. These involve the development, as the charged particles move away from the black hole, 
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of unspecified plasma-wave turbulence or strong internal shocks. Suggested shock generators are e + e~ pairs created 
by electromagnetic cascade in the postulated magnetic field of the emitted particles, or collisions of emitted particles 
with the interstellar medium. 

The Bclyanin, Kocharovsky, and Kocharovsky scenario makes the following assumptions: that the charged particles 
interact sufficiently to form a plasma; that a substantial magnetic field is generated by anisotropies in the charged 
particle distributions; and that significant turbulence arises in this field. In analyzing this scenario, we find that a 
number of the model components have been treated incorrectly. We conclude that a self-induced photosphere of this 
type should not occur around the black hole. 

First, Langmuir oscillations are longitudinal oscillations that arise in a plasma when the positively charged particles 
move collectively in the opposite direction to or independently of the negatively charged particles. The mutual 
attraction between the separated positive and negative charges provides a restoring force that constrains the charged 
particles to oscillate around their equilibrium positions in the rest frame of the plasma. The oscillations can propagate 
along this axis of particle motion. For example, Langmuir oscillations can occur in an ionic plasma placed in an external 
magnetic field. However, in the case of particles streaming away from a black hole, we note that there is no initial 
asymmetry in the distribution or response of positive and negative charges — for every particle, a neutral black hole 
stochastically emits on a short timescale an antiparticlc of equal mass and opposite electric charge 12 71 1 2811 . Any initial 
black hole charge is lost on a timescale much less than the age of the Universe for M^h < 1O 5 M [2-41. |25|. l2rl [27l. |28| . 

1/2 

The remaining random charge fluctuations are of order the Planck charge (he) = 11. 7e and therefore negligible 
[27L HH. Additionally, a neutral black hole has no intrinsic magnetic field. If emitted particles decay, total charge 
symmetry is preserved. In particular, the net electric charge of the final products of a QCD jet initiated by a decaying 
quark is equal to the negligibly small charge of the original quark, i.e. ±e/3 or ±2e/3. Moreover, self-induced QED 
turbulence is not observed in accelerator jets. Hence the outflow from a black hole remains essentially electrically 
neutral on all scales and Langmuir oscillations should not arise from charge asymmetries in the flux distribution. 

Although we do not expect a self-induced magnetic field to develop as the emitted particles stream away from the 
black hole, let us for the moment assume that collective wave-particle interactions and/or a turbulent magnetic field 
do arise in the charged particle flux when the hole is placed in an external magnetic field of sufficient strength and/or 
turbulence. Pions form via hadronization in QCD jets at a distance from the black hole of roughly 
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taking m q w 0.3 GeV for u and d quarks, with faster pions forming farther from the black hole and more coaxially 
than slower ones. The charged pions will decay at a distance 
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where E v = ^m^c 2 is the energy of the charged pion and t w ± = 2.6 x 10~ 8 s is the charged pion lifetime. The 
average energy of the pions produced by the black hole is [l2| 

E„ « 0.6\/WGeV GeV. (6.3) 

This value is distorted by the high energy E„ sa T particles but most pions in the instantaneous number flux are 
'slow' with energies £, k (1 - 3)m„-c 2 fl2l |. 

If Langmuir oscillations do arise in the charged pion flux, their frequency in the rest frame of the pion plasma is 
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where the prime indicates quantities in the plasma rest frame, e = 4.8 x 10 -10 esu is the electron charge in electrostatic 
units and n' n is the charged pion number density. Transforming into the rest frame of the black hole, the Langmuir 
frequency is 
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(see, for example, Rcf. |65] ) . This disagrees with Eq. (10) of Bclyanin, Kocharovsky, and Kocharovsky which omits 
a factor of 7^ 1 . In the black hole frame, the pion number density and flux are n T « 7V 7r ±/47rr 2 c and N^± « 
10 24 CWGcV) 372 s~\ respectively [3. 
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In the case of the slow pions with j w w 1 — 3, the first constraint for the hydrodynamical regime — that the Debye 
radius ~ c/u) n (which corresponds to r' DT . ~ c/w^, the distance at which the kinetic energy of a plasma particle 
balances its electrostatic potential energy in the plasma rest frame) be less than the characteristic size of the plasma 
r — is satisfied when Tbh ^ 1 GeV. For the fast pions with E n ~ E % , however, row is an order of magnitude greater 
than r at all Tbh and this constraint is never satisfied. In fact, at any given Tbh, the rjj^ < r constraint is satisfied 

1 /2 

only by pions with energy E n < 0.1 (Tbh/ GeV) ' GeV. Thus the claim by Belyanin, Kocharovsky, and Kocharovsky 
that rjjn < r is satisfied by all pions when Tbh is greater than a few GeV is not supported by our analysis. 

The second constraint for the hydrodynamical regime — that the Debye radius be greater than the average sepa- 
ration between charged pions in the plasma rest frame, i.e. t d-k > n ' n — corresponds to 

r Dv > 7 ^ 1/3 (6.6) 

in the black hole rest frame. The 7^ 3 factor is omitted by Belyanin, Kocharovsky, and Kocharovsky. In the case of 
the slow pions, this constraint is satisfied when the pions reach a distance from the black hole of 
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rw « ^ j cm. (6.7) 

In the case of the fast pions, it is satisfied when the pions reach a distance 

rw « 10- 15 f^M 1 cm. (6.8) 

Recall, however, that the fast pions never satisfy the rn ff < r constraint at any r. For all pions, we also have the third 
requirement — that the first and second constraints can be satisfied only if r is greater than the distance at which 
the pions form. Thus the region in which all three constraints are satisfied is max(r 7T i, r n7! ) < r < rv/ for slow pions 
and non-existent for fast pions. 

Repeating the above analysis for electrons and positrons, representative energies are E P O.iy/Tbh/ GcV GeV 
for the fast and E e « 130m e c 2 for the slow e ± [T§|- The total flux produced by a black hole is N e ± « 
10 24 (Ik/GcV) 372 s" 1 [13- For the slow e ± , it follows that r De < r when T bh > 10 3 GcV. For the fast e ± , 
TDe < t is never satisfied at any Tbh and the rjje ^5 T constraint is satisfied only by those e with energy E e < 
10 -3 (T(,h / GeV) 1 ^ 2 GeV. The second constraint, rjj e > r j^ 3 n e 1 ^ 3 1 is satisfied when the slow e ± have reached a 
distance r ne > 10~ 12 (Tbh/GcV) 1 ^ 4 cm from the black hole and when the fast e ± have reached a distance r ne > 

10~ 12 {Tbh I 'GcV) 3 ^ 2 cm. However, the e ± produced by the -k ± —> /i ± — > e ± decay chain form at a distance r e ; w 
10 5 7^ cm, which is much greater than r ne . Combining these results, the region in which all three constraints are 
satisfied is r e i < r for slow e ± . For fast e ± , there is no region in which all three constraints are satisfied. 

In the above discussion, we assumed that the strength of the ambient magnetic field was sufficient to generate 
hydrodynamical behavior. Let us now estimate the strength of the external magnetic field required to cause significant 
longitudinal separation of positive and negatives charges as they move radially away from the black hole. As we saw 
in Eq. (|4.6j) . the predominantly radially- moving particles may have a transverse velocity component vt ~ Tbhc/r in 
the black hole frame. Using the force equations F = ieurXfl and F = j 3 mdv/dt, the separation along the radial 
direction induced by the magnetic field is given by Ar^/r f=a eBrbh/') ' mc, implying 



Ars _ 1(T 26 / B \ ( Tbh \ f _ ^± 

r ~ 7 3 ^io- 5 Gy VGcvy 



for tt* (6.9) 



A,b_,10- 23 ^ B \fTbh_ fore± (610) 



r 7 3 V 10_5G / V GeV 

Certainly the Galactic magnetic field, which is of order 10~ 5 G, is too weak to produce significant charge separation 
and hence cannot generate Langmuir oscillations. Additionally, the magnetic field strength required for significant 
separation would appear to be unrealistically large in any astrophysical context. 

If we consider the pions, electrons and positrons to be generated in QCD jets of opening angle 0j et , the transverse 
velocity may be as large as vt ~ v6j e t ~ c9j e t and so the maximum separation is given by 



(6.11) 
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for tt^~ and 10 2 smaller for e^ . At accelerator energies, 9j et ~ 1.29a s [§6}. The maximum magnetic field-induced 
separation is then greater than the average separation between particles, n^ 1 / 3 , at 

10 7 9 / 4 / B y 3/4 (T bh \- 3 / 8 ± 

r ~^ v lb^Gj (.GevJ Cm *** (6 ' 12) 

jet 



> 



io 3 7 9/4 ( b y 3/i [T bh y 3 / 8 ± 

-^U^gJ UevJ Cm f ° re ' (6 ' 13) 



jet 



In this case, turbulence may be generated at high temperatures in a strong external magnetic field. However, only a 
small fraction of the emitted particles would have sufficient transverse velocity and the relevant orientation to satisfy 
the conditions to participate in Langmuir oscillations. We have already shown above that, once started, only the low 
energy ir^ and e ± can possibly sustain an MHD regime. 

We conclude that a self-induced magnetic field and MHD regime should not arise in the charged emitted particles 
as they stream away from the black hole. Hence no photosphere will be generated by this mechanism. If the black 
hole is embedded in an external magnetic field of sufficient strength and/or turbulence, we find that the MHD regime 
extends over a smaller region than Bclyanin, Kocharovsky, and Kocharovsky envisage and can involve only the lowest 
energy emitted particles, weakening the effect. Rees has also proposed a model in which the high energy electrons 
emitted by an expiring black hole form an expanding conducting shell which will be braked by any ambient Galactic 
magnetic field [67]. However, MacGibbon has found that the Rees mechanism, when re-analyzed including QCD 
decays, will produce a gamma-ray burst with a detection probability at Earth of much less than 1 (Tol. |64|. 



B. Cline and Hong Model 



D. Cline, Matthey and Otwinowski have proposed that Galactic PBHs may explain the shortest-lived gamma-ray 
bursts Q. They suggest two mechanisms for burst generation: plasma formation at T^k ~ 10 GeV — 10 TeV, together 
with the MHD instability proposed by Belyanin, Kocharovsky, and Kocharovsky [22]; or 'fireball' formation at the 
QCD threshold Tbh ~ Aqcd, as proposed by Cline and Hong [1,0, HI]. As we have seen, the Belyanin, Kocharovsky 
and Kocharovsky mechanism is not viable in the absence of a large external magnetic field. The second mechanism 
postulates that the emitted particles undergo a first-order phase transition around Tbh ~ Aqcd, accompanied by rapid 
loss of the remaining black hole mass via a Hagedorn-type exponential increase in the number of hadronic degrees of 
freedom at Tbh ~ Aqcd, without the black hole temperature climbing above Kqcd- These hadrons are assumed to 
be directly emitted by the black hole, and not the result of decays. 

The Hagedorn model was suggested by the exponential increase in hadronic resonances observed in early accelerator 
experiments at these energies. However, the resonances have long since been understood in terms of the more 
fundamental quark model, and accelerator energies well above Aqcd have long since been achieved. The standard 
Hawking picture, in which the black hole emits only those particles which are fundamental on the scale of the black 
hole and which was discussed in detail in Section IV. K, is consistent with current accelerator observations. Thus a 
Tbh ^ Aqcd hole should directly emit quarks and gluons, which then hadronize at a distance r w Tbh/ Aq CD from 
the black hole, rather than an explosive Hagedorn-type particle distribution. Similarly an m T < Tbh ;$ Aqcd black 
hole should emit pions as fundamental particles. Furthermore, in the Hagedorn picture the exponential distribution of 
resonances are possible particle states but whether the states are occupied is determined by the available energy. As 
we have seen in Section IV. K, the energy per emitted quark or gluon around Tbh ~ Aqcd is only E ~ Aqcd- Dense 
occupation of the Hagedorn resonances would require a direct coupling of the black hole mass to the Hagedorn states. 
If the black hole is rapidly losing mass into the Hagedorn spectrum, however, Eq. (|2.2p dictates that the temperature 
must rapidly climb above Age d ■ This would violate the condition of the original Hagedorn model [69j in which the 
exponential growth in the density of states results from the temperature being held constant. 

Hence, based on the observed behavior in accelerators, there should be no rapid increase in the black hole evaporation 
rate as Tbh — * Aqcd and no gamma-ray photosphere burst of the type envisioned by Cline and Hong. The 'fireball' 
model proposed by Moss [7(j , in which a phase-transition at the Aqc d threshold generates a confining bubble around 
the black hole, is also not consistent with accelerator observations. Note that the Cline and Hong 'fireball' model is not 
related to the Heckler brcmsstrahlung-induccd photosphere. A Heckler-type photosphere involves only interactions of 
particles after emission and is not accompanied by a rapid increase in the Hawking evaporation rate. 

Although our analysis rules out photospheres in such models, we cannot exclude the possibility that black hole ex- 
plosions could generate short-period gamma-ray bursts by some other as yet unconsidered mechanism or by placement, 
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for example, in unusual magnetic field configurations. The possibility that one might detect black hole explosions (and 
thereby test the Hawking theory) at the present epoch is so important that all such scenarios are worth exploring. In 
particular, it is important to examine whether short-period gamma-ray bursts have other characteristics expected of 
a PBH population, independent of knowledge of the specific burst production mechanism. 

VII. CONCLUSIONS 

After careful examination of the Heckler photosphere model, we find that the emitted particles are not expected 
to interact sufficiently to form a QED or QCD photosphere around an evaporating black hole. Key to determining 
the fate of the emission is understanding the causality constraint and scatterings which occur within a distance 
~ l/m e or rsj 1/Aqcd of the hole. Our main point is that no QED photosphere can form around a black hole 
because a particle must be causally connected to any particle it interacts with and requires a distance ~ E/m^ to 
form in each bremsstrahlung interaction. This means that each particle is highly unlikely to interact and generally 
undergoes at most one net interaction as it streams away from the hole. Neither Heckler nor Clinc, Mostoslavsky, and 
Servant include these effects. Instead they envision a scenario in which each emitted particle undergoes a sequence of 
interactions. In Heckler's model these start at a radius less than rv em and continue to a radius greater than rv e mj m 
the Cline, Mostoslavsky, and Servant simulations, they start at a radius much less than r^ em and continue to 7V em . 
Furthermore, the modification required to their models is not a simple one, since the causality and formation distance 
constraints and a corrected geometrical description must go into the analysis at the start and cannot simply be added 
later. Additional QCD phcnomenological arguments, based on the limited energy per initial particle and availability 
of QCD final states, also rule out QCD photosphere (chromosphere) development around Tbh ~ h-QCD- 

Although we find that photosphere formation is not supported, an emitted particle may occasionally undergo one 
net QED interaction within a distance ~ l/m e of the black hole or one net QCD interaction within ~ 1/Aqcd. 
However, this interaction could only very slightly modify the observable 7, p, p, e + , e~, v and v spectra from a 
PBH background or an individual high-temperature black hole. The spectra remain essentially those of the standard 
Hawking model, with little change to the detection probability. We also find that interactions are insufficient to 
support other models for photosphere formation or enhanced PBH signals proposed in the literature, unless the black 
hole is located in an unusual ambient environment or Standard Model physics breaks down. 

In this paper, we have not addressed the emission expected from higher-dimensional black holes which may be 
created at TeV or higher energies in accelerators [zj 0, [TH, [71, IzHl or cosmic ray interactions [7(| [73] ■ Our analysis 
implies that the spectra which were calculated, in for example Refs. [78|, [79| . assuming Heckler- type interactions 
between the particles emitted in the final burst of evaporation from a higher-dimensional black hole will need to be 
modified. In particular, the issues raised in Sections IV. B, D, F, G, H and K will need to be applied to the analysis 
of the spectra from higher-dimensional black holes. In the case of accelerator experiments, the relevant condition 
to investigate is whether M ~ 1 is approached, because even one interaction per particle may modify the spectra 
significantly compared with experimental precision. A number of the issues we raise were independently recently 
included in the numerical simulation of the decay of CERN LHC TeV-scale black holes by Alig, Drees and Oda 
[80l ] . They concluded that scattering between partons emitted by the higher-dimensional black hole is insufficient to 
form a thcrmalizcd chromosphere. Although their methodology was restricted to scattering between partons which 
approach each other and included a lower cutoff of 1 GeV in the exchanged transverse momentum, their results provide 
insight into the applicability of the Heckler scenario to TeV-scale higher-dimensional black holes. Our discussion of 
geometrical and causality considerations in Section IV. F and in our accompanying paper [39j may be viewed as a 
partial justification of the approximation used in [8(| of only including scattering between partons which approach 
each other. 
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